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In this paper we reconsider the model of neutrino production during the 'bright phase', first 
suggested in 1977, in the light of modern understanding of both the role of Pop III stars and of 
acceleration of particles in supernova shocks. We concentrate on the production of cosmogenic ultra- 
high energy neutrinos in supernova (SN) explosions that accompany the death of massive Population 
III stars. Protons are assumed to be accelerated at such SN shocks and produce neutrinos in 
collisions with CMB photons. In the calculations we deliberately use simplified assumptions which 
make the physical results transparent. Pop III stars, either directly or through their SN explosions, 
are assumed to be responsible for the reionization of the universe as observed by WMAP. Since 
the evolution of massive Pop III stars is much faster than the Hubble rate H , we consider the 
burst of UHE proton production to occur at fixed redshift [zb = 10 and Zb = 20). We discuss in 
some detail the problems involved in the formation of collisionless shocks in the early universe as 
well as the acceleration of charged particles in a medium that has potentially very low pre-existing 
magnetization, if any at all. The composition of the accelerated particles in Pop III stars explosions 
is expected to be proton dominated, based upon the predictions of BBN and the Hydrogen-enhanced 
stellar-wind from primary Pop III stars. A simple calculation is presented to illustrate the fact that 
the diffuse neutrinos flux from the bright phase burst is concentrated in a relatively narrow energy 
interval, centered at El = 7.5 x W 15 (20/ Zb) 2 eV. The flux may be detectable by IceCube without 
violating the cascade upper limit and without exceeding the expected energetics of SNe associated 
with Pop III stars. A possible signature of the neutrino production from Pop III stars may be the 
detection of resonant neutrino events (i> E + e~ — > W~ hadrons) at energy _Bo = 6.3 x 10 eV. 
For the burst at Zb = 20 and ^ e -flux at the cascade upper limit, the number of resonant events in 
IceCube may be as high as 10 events in 5 years of observations. These events have equal energies, 
E = 6.3 x 10 15 eV , in the form of e-m cascades. Taking into account the large uncertainties in the 
existing predictions of cosmogenic neutrino fluxes at E > 10 15 eV, we argue that UHE neutrinos 
from the first stars might become one of the most reliable hopes for UHE neutrino astronomy. 



I. INTRODUCTION 

Ultra High Energy (UHE) neutrino astronomy at ener- 
gies above 10 16 — 10 17 eV may open a new observational 
window on the universe. The most reliable prediction for 
the diffuse fluxes of these neutrinos can be achieved for 
the cosmogenic neutrinos [l|, produced inpy collisions 
of UHE protons with CMB photons (see || and 3 for 
recent calculations and reviews). Such reliability mainly 
follows from three factors: i) The flux of parent protons 
at redshift z of neutrino production can be normalized to 
the observed flux of protons at present (redshift z = 0); 
ii) The number density and energy spectrum of the tar- 
get (CMB) photons are both known with high accuracy 
at any cosmological epoch z; Hi) the interactions of UHE 
protons with CMB photons occur at center of mass (cm) 
energy that does not exceed ~ GeV, well accessible to 
laboratory experiments and hence well known. 

On the other hand, the numerical predictions of neu- 
trino fluxes at Earth are affected by huge uncertainties, 
mainly due to the unknown evolution with redshift of 
the source space density and luminosity, the generation 
spectrum (oc E~ lg ) and the maximum energy (E^ x ) of 
accelerated protons. Because of these uncertainties, any 
upper bound on the flux of cosmogenic neutrinos becomes 
of the highest importance. For cosmogenic neutrinos this 
limit is not provided by the Waxman-Bahcall bound [|J, 



because the latter is imposed by the parent proton flux 
at z = 0, fixed in the calculations. Meanwhile the cal- 
culated flux can differ by orders of magnitude due to 
uncertainties in 7 g , E ma _ x and the evolution of sources 

a. 

A more general upper limit, valid in particular for the 
cosmogenic neutrino flux is the cascade upper limit 
The production of UHE neutrinos from pion decays is ac- 
companied by the production of UHE electrons, positrons 
and photons, which start an e-m cascade due to the colli- 
sions with CMB photons. The HE photons, remnants of 
these cascades, contribute to the diffuse gamma-radiation 
with gamma rays of energy reaching several hundreds 
GeV. Recent measurements of the diffuse gamma-rayflux 
by the Fermi-LAT telescope Q was used in Ref. [H to 
put a stronger upper limit on the diffuse neutrino flux, 
which excluded many models with detectable cosmogenic 
neutrino flux. 

The detectability of cosmogenic neutrinos is fur- 
ther disfavoured by the recent Auger measurements of 
UHECR mass composition 0] : both X max and RMS data 
strongly suggest that at energies larger than 2 x 10 18 eV 
the primaries are nuclei with steadily increasing mass 
number A. In contrast with these data the HiRes [lfj and 
Telescope Array [ll| mass composition agrees well with a 
pure proton composition. The cosmogenic neutrino flux 
produced by primary nuclei is considerably lower than in 
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the case of primary protons [12j. The Auger results on 
mass composition, together with the Fermi-based neu- 
trino upper limit suggest a rather grim possibility of de- 
tection of cosmogenic neutrinos. 

Here we revise a model of UHE neutrino production 
related to the so-called Population III (Pop III) stars 
[T3I [LH ] . Recently this population of stars attracted much 
attention as the most plausible sources of ionizing pho- 
tons in the early universe, at redshift z ~ 10 — 20, as 
observed by WMAP 

Originally, one of the main motivations that led to pro- 
pose the existence of Pop III stars was the apparent gap 
between the chemical composition of the universe as pre- 
dicted by BBN, mainly light elements up to Lithium, 
and the presence of metals needed for the formation of 
Population II stars. Pop III stars bridge this gap by syn- 
thesizing the first metals in the universe and spreading 
them through the interstellar and intergalactic medium, 
where they are later instrumental for the formation of 
Pop II stars. For this scenario all Pop III stars or a large 
fraction of them must finish their evolution by SN explo- 
sions, and therefore these stars should be very massive, 
heavier than 50M Q . 

After recombination at z rcc rs 1100, the baryons in the 
universe are in the form of cold, neutral, atomic hydrogen 
and light atoms. However, at redshift z rc ion = 11.0 ± 1.4, 
according to WMAP [l|, the baryonic gas is observed 
to be ionized. The reionization occurred due to photons 
produced by astrophysical sources. The most plausible 
sources of such photons are indeed hot Pop III stars, 
which produce ionizing photons either directly [H| or as 
a result of Pop III supernova explosions [TtJ • 

Pop III stars arise from baryons gravitationally pulled 
into potential wells formed by dark matter (DM) dom- 
inated mini-halos. These mini-halos collapse at z ~ 
20 — 30 and eventually capture baryons which can form 
stellar objects. Numerical simulations of this collapse 
show no fragmentation of baryonic matter and the first 
stars form with large masses, spread over a wide range, 
with a typical value M ~ 100M Q (see [l8| and references 
therein) . 

At masses 140M o < M* < 260M©, Pop III stars 
are fully disrupted in SN explosion due to e + e~ pair- 
instability. These SNe are called Pair-Instability SNe 
(PISN). The energy output of PISN in the form of ejecta 
is in the range 10 51 — 10 53 erg, and could be much higher 
than for Pop II and Pop I SNe. In the mass range 
M < 140M© and M > 260M Q , a massive black hole 
and disc are produced, with ejection of gas in the form 
of jets. These stars are good candidates for GRBs, and 
UHE neutrinos from Pop III GRBs were recently studied 
in [19[ . For further information on Pop III stars and SNe 
see the reviews in [20j |. 

In the following we will refer to this epoch of Pop III 
star formation as the bright phase, an expression often 
used also to refer to a period of enhanced quasar activity 
(about ~ 20 quasars are observed at 5.7 < z < 6.4 (2l|). 

The paper is organized as follows: in section |H] we dis- 



cuss some subtle issues related to the formation of shocks 
in Pop III supernova explosions and in the particle accel- 
eration at such shocks. In section Mil the basic features 
of the bright-phase model are described. In section IIVI 
we present the calculations of the neutrino fluxes and the 
constraints on such fluxes. The conclusions are presented 
in section IVl 



II. SHOCK FORMATION AND PARTICLE 
ACCELERATION 

The environment in which Pop III stars explode as SNe 
is quite unlike that of SNe in the present universe and one 
may wonder whether particle acceleration may take place 
at all. The intergalactic medium at redshift z ~ 10 — 20 is 
unlikely to be magnetized in any appreciable way. This 
leads to two questions: 1) does a shock develop in the 
supernova explosions, since such shocks are collisionless, 
namely mediated by electromagnetic instabilities? 2) if 
the shock forms, is there enough magnetic turbulence 
to lead to diffusive particle acceleration up to energies 
~ 10 19 — 10 20 eV, needed for neutrino production? 

Numerical simulations of the development of the ion 
Weibel instability in the absence of a pre-existing mag- 
netic field show that a shock front does in fact form [22| • 
The simulations are carried out for relativistic electron- 
ion plasmas, while the motion of the plasma ejected from 
a Pop III star explosion might be newtonian or relativistic 
depending upon details of the explosion, but the general 
physical principles that lead to the formation of the shock 
front should be left unchanged, so we may be confident 
that a shock wave does indeed form even in a weakly 
magnetized or unmagnetized medium such as the one 
that is expected to embed Pop III stars. The fate of 
the magnetic field behind the shock is not well estab- 
lished: the magnetic filaments which are formed through 
the Weibel instability might merge to form larger scale 
magnetic fields, but it is not clear whether this happens 
before the field may be damped. 

More important than the presence of pre-existing mag- 
netic field is the reionization of the surrounding gas prior 
to the SN explosion: a collisionless shock would hardly 
develop in a tenuous neutral medium, such as in the uni- 
verse soon after recombination. Pop III stars can natu- 
rally fulfill the ionization requirements since, being very 
hot, they are expected to be responsible for the produc- 
tion of intense UV radiation that causes reionization (see 
e.g.[16| and references therein). 

The issue of particle acceleration is more delicate in 
that diffusive acceleration at a shock front requires the 
existence of a turbulent magnetic field upstream and 
downstream of the shock. The only way that magnetic 
field can be generated upstream of the shock is through 
streaming instability, an intrinsically non-linear process 
(particle acceleration occurs if accelerated particles are 
able to excite the instability but the instability is excited 
by accelerated particles). This process is also necessary 
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in SNe in the ISM, and in fact the instability must be 
driven to its non-linear regime if to reach cosmic ray en- 
ergies as high as the knee 12311. The streaming instability 
can proceed in a resonant |24{ or non-resonant (25j way, 
but these are just two sides of the same physical phe- 
nomenon [2(| . The non-resonant branch grows faster for 
higher shock velocities [261 ] . 

Once the instability turns non-linear it is very difficult 
to predict its saturation, which might be dominated by 
intrinsic dynamical scales (such as the advection time of 
a fluid clement through the shock) or by damping. The 
issue of damping might be of particular importance in 
the clouds where Pop III stars originate since the neu- 
tral fraction might be appreciable and ion-neutral damp- 
ing could well shut the instability off. In the absence of 
damping and assuming a naive extrapolation of quasi- 
linear theory to the non-linear regime one can estimate 
the saturation levels of the instability as: 
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8tt M a 
in the resonant case, and 



4tt 2 P s c^ CR 



(1) 



(2) 



in the non resonant case. Here Ma = V s /v a is the 
Alfvenic Mach number, and v a = B /(iirp) 1 ^ 2 is the 
Alfven speed. Clearly neither the Alfvenic Mach number 
nor the Alfven speed are well defined in the non linear 
regime, but this ambiguity is the price to pay to write 
down a simple estimate of the magnetic field generated 
through streaming of accelerated particles upstream of 
the shock. In all these formulae we assumed that the 
shock is non-relativistic. Our estimate for the amplified 
magnetic field is: 



6B res K96B^n\ / X / % / R»G, 
in the resonant case and 

6B nr a 590n^V 9 3/2 ^ /iG, 



(3) 



(4) 



in the non resonant case. Here Vg = Vs/10 9 cm/s and 
n\ = n/lcm~ 3 . ^cb, is the fraction of the inflow ram 
pressure m p nV 2 that gets converted to accelerated par- 
ticles. 

Assuming Bohm diffusion, the acceleration time in the 
two cases reads 



(5) 



r res (E) « 3.5 x 10- 7 E(eV)B-y 2 n- l/4 V 9 - 5/2 ^R / - 

in the resonant case, and 

r nr (E) » 5 x 10- s E(eV)Vi 7/2 n; 1/2 ^ 2 s (6) 

in the non-resonant case. 

It is important to keep in mind that in the non- 
resonant case the self-generated turbulence is on spatial 



scales that are much smaller than the gyration radius of 
the particles, therefore these modes are not very effective 
in scattering the particles, despite the large growth rate. 
In this sense, the assumption of Bohm diffusion for non- 
resonant modes is poorly justified, and would require an 
efficient cascade of the modes to much larger scales, an 
inverse cascade. 

There are at least three different ways to define the 
maximum energy of accelerated particles: 1) Time limi- 
tation: t(E^ x ) ~ Tst, where T$t is the time when the 
remnant enters the Sedov- Taylor phase; 2) Space limita- 
tion of the precursor: D{Ef^ x )/V s ~ tjRst, where Rst 
is the radius of the supernova remnant at the beginning 
of the Sedov phase, and r\ < 1; 3) Space limitation on 
the gyration radius: rtCEwJ ~ Rst- The three crite- 
ria lead to different estimates for the value of E^ x . 

The Sedov time and radius are related to the ejected 
mass and to the total energetics through: 



M. 



5/6 
ej',0 



and 



T ST = 70 i/ 2 T/3 
e 51 n l 



Rst = 6.6 x 10 18 ( 



(7) 



(8) 



where €51 is the total energetics in units of 10 erg and 
M e j,Q is the mass of the ejecta in units of solar masses. 

Among the three criteria for -E^ax listed above, the 
last one leads to the highest estimated value: 



E acc ^2xW 17 B 1 / 2 n- 1/12 V 1/2 £ 1/2 M 1/3 eV 



for the resonant case, and 



pa 



10 1 



S 1/6^3/2^1/2 1/3 
n l V 9 t iCR 1V1 ej,e 



eV 



(9) 



(10) 



for the non resonant case. One can see that maximum 
energies as high as (2 — 5) x 10 19 eV might be reached if 
the values of the parameters are pushed to their extremes. 

In principle larger values of -E^ax can ^> e reached if 
the Pop III stars explode in a dense stellar region, so 
that accelerated particles may feel the repeated action of 
multiple shocks. Such situation occurs at high redshifts 
where the physical volume which hosts the fixed number 
N of pregalactic Pop III stars is (1 + z) 3 times smaller 
than at present. Even a stronger case might be realized in 
mini-cluster models, when the Pop III stars are produced 
in the relatively small volume of a mini-cluster. In this 
case the spatial scale to be used to determine the max- 
imum energy is the size of the region where the shocks 
emerge, provided the acceleration time remains smaller 
than the loss time scale as plotted in Fig. [TJ Acceleration 
at multiple shocks has been first discussed in |24| and is 
known to result in a hard spectrum: in the asymptotic 
case of a very large number of shocks, the spectrum of 
accelerated particles tends to cx E~ 3 / 2 . 
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Larger values of Ef^, up to 10 20 eV, were inferred in 
previous literature for the case of acceleration at relativis- 
tic shocks with large Lorentz factors T >• 1, see e.g. [13] ■ 
Usually this case is considered for jets, and in particu- 
lar for GRBs, which are often assumed to be potential 
sources of particles with maximum energies ~ 10 20 eV 
[IM [29| . These energies are the consequence of the rela- 
tivistic motion of the outflow in jets, which in principle 
can shorten the acceleration time by a factor 1/r. 

Particle acceleration at relativistic shocks requires the 
presence of strong turbulence generated in the down- 
stream region, in order to avoid particle trapping, typi- 
cal of relativistic shocks: particle acceleration is possible 
only for quasi-parallel shocks namely for magnetic fields 
oriented within an angle ~ 1/r from the normal to the 
shock surface. At perpendicular shocks the return prob- 
ability from downstream tends to vanish thereby making 
spectra steeper. In fact, even the compression of the large 
scales in the upstream turbulent fields at the shock sur- 
face may lead to spectra steeper than the canonical E~ 2 3 
[30| , that is expected for parallel relativistic shocks in the 
regime of small pitch angle scattering. The presence of 
strong turbulence may alleviate this problem, so to make 
GRBs from Pop III stars potential sources of UHECRs 
with energy £^ c a c x ~ 10 20 eV. 

The issue of whether particle acceleration can occur 
in the first stars is also of relevance for the origin of 
cosmic magnetic fields: in a recent paper [3l[ the au- 
thors proposed that if supernovae arising from the death 
of primeval stars accelerate CRs effectively (though not 
necessarily to very high energies) , the instability induced 
by the escape of these cosmic rays into the intergalac- 
tic medium may lead to the formation of magnetic seeds 
that can possibly be reprocessed and amplified at later 
cosmic epochs. 



III. THE MODEL 

Here we consider a model in which the bright phase is 
powered by SN explosions of Pop III stars. 

In order to enrich space with metals, the rate of SN ex- 
plosions must be high and thus Pop III stars have to be 
massive, typically with M > 50M Q . Stars in the interval 
of masses (140 — 260)M Q undergo e + e~ instability and 
end their evolution with full destruction in SN explosions. 
Outside this interval, and most notably at M > 50M Q , a 
SN explosion leaves behind a massive black hole. In both 
cases a SN explosion results in the formation of a shock 
front where particles from the interstellar medium can 
be accelerated (see section [TT] for the detailed discussion). 
The formation of such (collisionless) shock is expected to 
be made possible because of the ionization of the circum- 
stellar region caused by photons produced by either the 
Pop III star itself or its supernova. The chemical com- 
position of the accelerated particles should reflect that of 
the circumstellar medium, which is quite different from 
that in our Galaxy. The primordial gas after BBN is 
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FIG. 1: Energy losses of protons at cosmological epochs 
Zb — 10 (shown in red) and Zb ~ 20 (shown in blue). Adia- 
batic energy losses are given by H(z). Also shown are losses 
due to pair production on CMB, p + 7 cm b — > p + e~ + e + , 
and pion photo-production, p + 7 cm b - ► iV+pions. Photo- 
production responsible for the generation of neutrinos dom- 
inates at energy above the crossing of the pair-production 
and pion production curves (E b c = 3 x 10 18 eV at z b = 20 
and E\ = 5.5 x 10 18 eV at z b = 10). The intersection of 
the pair-production curve with H{z) (adiabatic losses) de- 
termines the energy e pa ir, above which the pair-production 
energy losses dominate (e pa ir = 2.1 X 10 16 eV at z b = 20 and 
Epair = 5 x 10 16 eV at z b = 10). 



characterised by 75% of Hydrogen and by absence of the 
heavy elements (metals). Some heavy element pollution 
is expected to be caused by multiple SN explosions in the 
same spatial regions. On the other hand, the metallicity 
of the gas in Pop II star formation regions is expected 
to be relatively low, therefore it appears reasonable to 
assume that the medium around Pop III stars is domi- 
nated by light elements. We assume here that most of 
the accelerated particles are protons. 

As discussed above, the SN explosion of Pop III stars 
can result in the formation of either non-relativistic or 
relativistic shock fronts, the latter usually associated to 
the formation of jets. The spectrum of accelerated parti- 
cles in the two cases is usually taken to be cx E~~< 9 with 
7 g = 2.0 for non-relativistic shocks and 7 g = 2.2 — 2.3 
for relativistic shocks, although numerous physical effects 
can change this naive expectation. In both cases we as- 
sume maximum energy of acceleration E^ x ~ 10 20 eV 
and the minimum energy to be as low as E min ~ m p c 2 ~ 
10 9 eV. Increasing E m i n results in larger neutrino flux in 
the calculations below. 

As far as E max is concerned, the energy 10 20 eV may 
be problematic in non-relativistic shocks, although par- 
ticle reacceleration at multiple SN shocks may hopefully 
provide this energy (see section [II]) . The situation with 
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E max ~ 10 20 eV is more promising but still poorly un- 
derstood in the case of relativistic shocks. 

Rather than achieving a quantitatively accurate pre- 
diction of the neutrino flux from Pop III stars, our aim 
here is to develop a simple calculation that catches the 
main physical ingredients and illustrates in a clear, trans- 
parent way, the relevant aspects of the problem. The dis- 
tinctive feature of our model is the assumption that the 
duration of the bright phase, occurring at redshift z b , 
is very short compared with the Hubble time H^ 1 {z b ) 
at the same e poch. The latter can be calculated as 
H{z b ) = H oy /n m (l + z b ) 3 + Q A with H = 72 km/s 
Mpc, fi m = 0.27, and fi A = 0.73. At z b = 10 
the Hubble time is H(z b )^ 1 ~ 8 x 10 8 yr, to be com- 
pared with the lifetime ~ 10 6 yr of a Pop III star with 
M ~ IOOMq on the Main Sequence (the longest time 
of this star evolution) , with the time of particle accelera- 
tion, t acc ~ 10 4 — 10 5 yr, and with the characteristic time 
of photopion energy losses, ~ 3 x 10 5 yr, responsible for 
neutrino production (see Fig. [1] for comparison of these 
times at various energies). These estimates serve as a 
clear justification for the assumption of bursting bright 
phase: all processes are assumed to occur at the fixed 
redshift z b , for which further on we consider z b = 10 
and z b — 20 as our benchmark cases (see Fig. Qp. It is 
interesting to note that models with two bursts of reion- 
ization can be formulated |32j, as based on taking into 
account the important role of neutral H2 molecules in 
star formation. The redshifts of the predicted bursts are 
Z\ ~ 11 — 15 and z^ ~ 6, close to the value of z b that we 
consider. 

Having in mind the same goal of simplicity of formulae 
and transparency of the calculations, we assume power 
law spectra with slope 7 9 = 2.0 , which does not lead to 
dramatic differences compared with similar spectra with 
7 ff ~2.2-2.3. 



IV. ASSUMPTIONS AND ESTIMATES 

The burst of accelerated protons produced during the 
bright phase at epoch z b produces a spectrum that can 
be expressed as a function of the energy density oj p (z b ) 
of the same protons at that epoch as: 



n p (E b ,z b ) 



uj p (z b ) 



i„ ( 77* max / rp 

tn(h b /E t 



mm \ 



-e: 



(ii) 



Henceforth uj p (z b ) will be used at the main normalization 
parameter of our calculations. 

The cosmogenic neutrinos are produced in collisions 
of UHE protons with CMB photons through generation 
of pions, p + 7cmb — > 7r ± +anything. On average a pro- 
ton transfers 20% of its energy to the leading pion, and 
this energy is approximately equally divided among the 
four leptons (three neutrinos and one electron), so that 
E u = (1/20)£7 P . Photopion production dominates start- 
ing at the proton energy E^, where energy losses due 
to photopion production exceed pair-production losses 



(p + 7cmb — > p + e~ +e + ). At z = this condition occurs 
at proton energy E c » 6 x 10 19 eV, therefore at epoch z b 
this energy is E^(z b ) = E c j (1 + z b ) (see Fig. [J), since the 
energy of CMB photon is (l + z b ) times higher. Hence the 
characteristic neutrino energy in the observed spectrum 
is 



El 



1 



E r . 



20 (1 + z b f 



7.5 x 10 1 



20 



1 



eV, (12) 



where an extra factor (1 + z b ) appears because of the 
redshifted neutrino energy. Below E°, the neutrino spec- 
trum becomes flatter, being produced by pions moving 
backwards in the cm-system. Above the neutrino 
spectrum follows the parent proton spectrum, cx E~ 2 . 
Eq. (fT2|) thus gives the neutrino energy where detectabil- 
ity reaches its maximum. The main contribution to the 
number of detected neutrinos is given by some energy 
interval centered at E^. 

For the calculation of the neutrino spectrum it is con- 
venient, following [3, [33, Hf| , to introduce the so-called 
'unmodified' proton spectrum, calculated taking into ac- 
count only adiabatic energy losses due to redshift. This 
proton spectrum at z = is 



\E) = 



io p (z b ) 



1 



(1 + z b ) 4 \n(E max /E n 



-E~ 



(13) 



Eq. (fT3|) is easy to understand because uj p (z b )/(l + z b ) 4 
is the unmodified energy density of protons at z = 0: 
(1 + z b ) 3 is due to the expansion of the universe and the 
extra factor (1 + z b ) accounts for adiabatic energy losses. 

The spatial density of UHE neutrinos at z = can 
now be readily calculated using the unmodified proton 
spectrum and taking into account that each proton in 
one photonuclear pj collision produces three neutrinos 
and the energy of a neutrino E v is connected with the 
energy E of the 'unmodified' proton as E v = (1/20).E, 
independent of the redshift z when the collision occurs, 
because adiabatic energy losses of neutrinos and UHE 
protons are identical. 

Thus, one obtains at E v > El : 



n v (E v )dE v 



3/ ff ± uj p (z b ) 



1 dE, 



20 (1 + ztf MEmuc/EmSn) 



r. ( 14 ) 



where f n ± w 2/3 is the fraction of charged pions. 

The maximum neutrino energy at z = can be esti- 
mated as 



1 



1 E^ 



20 max 201 + z b 



The neutrino flux, 



E?,J, 



4tt (1 + z b y HE max /E min ) ' 



(15) 



(16) 



is fully determined by the value of the basic parameter 
oj p (z b )/(l + z b ) 4 . 
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We can normalize this parameter to the future IceCube 
sensitivity [H E 2 J U ^ (E)=3x 1(T 9 GeV/cm 2 s sr. (The 
present upper limit of IceCube is a factor of 3 higher). 
From this condition the basic parameter must be 



U p {Zb)/(\ 



z b f > 9.5 x 10~ 7 cV cm" 



(17) 



We can estimate now the maximum energy of acceler- 
ation = _E" lax discussed in section [TTJ The rigor- 
ous limit follows from the condition that pion production 
dominates upon pair production (see Fig. [T]): 



1 



Zb 



3x 10 



IS 



20 



1 



Zb 



eV. (18) 



The same limit follows from the condition which provides 
an appreciable neutrino flux, _E™ ax > E%. In the present 
estimates we assume that £!™ ax ~ 10 20 eV. The case 
E™ 3 *- > 1 x 10 19 eV needs more detailed calculation of the 
neutrino spectra and will be presented in a forthcoming 
paper. 

In the following we discuss the constraints on the cos- 
mogenic neutrino flux from the bright phase. 

The production of UHE neutrinos from pion decays is 
accompanied by e-m cascade radiation, which provides 
us with an upper limit on the neutrino flux [6j. The 
neutrino flux given by Eq. (|16[) with energy density (|17[) 
must respect this limit as given by the recent Fermi ob- 
servations • In @ this limit is expressed in terms of an 
upper limit on the energy density of the cascade radiation 
w™| x <5.8xl0~ 7 eV/cm 3 . We can estimate the cascade 
energy density in our calculations from u) p {zb) evaluated 
in Eq. (17]). 

At energies Eb > e pa ir (see FigHJ, protons lose energy 
on time scales shorter than the Hubble time H^ 1 (z b ), 
producing e + e~ pairs and pions. The energy density w£ as 
of the cascade radiation at epoch Zb can be written as 



cas 



n p (E b )E b dE b 



(19) 



and for the cascade energy density at z = we have 

-l 
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Up(z b ) 
(1 + z b ) 4 
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^max 



= 2.9 x 10 _7 eV/ 



7^111111 
^b 

-7„ 



(20) 



One can see that the cascade energy density (120)) . which 
follows from neutrino production in the bright phase does 
not exceed the upper limit w™ ax = 5.8 x 10~ 7 eV/cm 3 , 
obtained in [8j from Fermi observations. In fact, this up- 
per limit is derived for the cascade spectrum produced at 
redshift z — 0, while in our case the cascade is produced 
at large Zb and its spectrum is redshifted. This makes 
the upper limit higher than w™ ax from Q. 

Finally, we address the issue of whether the energy 
density of protons in Eq. (fl"7| which is needed to warrant 
detectability of the neutrino flux in IceCube can be pro- 
vided by Pop III SN explosions. For this, we estimate a 



fraction £ of cosmological baryon mass processed by Pop 
III pre-supernovae. 

From the energy density of UHE protons u> p (zb) — 
Wg N n^(zb), provided by the SN energy release in the 
form of UHE protons Wg N and from the space density 
of Pop III presupernovae n* , one obtains 



Wp(z b ) M*(l + Zb) 

(i + zbY wg N n bPci ' 



(21) 



Here we used the equality M*n*(zb) = £pb(zb), where 
A/* ~ lOOAf© is a typical mass of Pop III presupernova 
and pb(zb) = f2fc/9 C r(l + Zfc) 3 is the total baryonic mass 
density at epoch Zb, where for the baryonic gas density 
at z = we use the WMAP pi value 6 p CI = 4.27 x 
10~ 31 g/cm 3 . For M» = 100M© the SN energy release 
can be conservatively estimated as Wg N ~ 10 51 erg. As 
a result we obtain for the case of Eq. (|T7|) the reasonable 
value £ = 1.4 X 10~ 2 , i.e. about 1% of the total baryonic 
mass processed by Pop III presupernovae. 

If to assume that each SN explosion leaves behind 
a black hole with mass Mbh ~ 3M e , the fraction of 
baryonic mass in the form of these black holes is only 
rj = £M bh /M* - 4 x 10- 4 . 



Ve + e 



hadrons resonance. 



The resonant production of hadrons by UHE v e - 
neutrinos is a remarkable feature and signature of the 
Pop III bright phase model. 

The resonant production of W-bosons was first pro- 
posed by S.L. Glashow in 1960 [3(| for the generation of 
high energy muons in the reaction v e + e~ — > W~ — ¥ 
fi~ + Va- In 1977 Berezinsky and Gazizov (3?| suggested 
another channel of the Glashow resonance v e + e~ — > 
IU~ — > hadrons for the detection of UHE neutrinos. The 
resonant neutrino energy is given by 



Eo = m 2 w /{2m, e ) = 6.3 x 10 6 GeV, 



(22) 



and the rate of resonant events ^ rcs in a detector with 
total number of electrons N e is determined by the exact 
formula derived in 1371 : 



2tt N e a cff J Pc (E )E 



(23) 



where Jv k [Eq) is the diffuse flux at the resonant energy 
Eq, N e — (5/9)M/mu is the number of electrons in the 
underground detector with mass M, 2tt is the solid angle 
at which a deep-underground detector is open for res- 
onant neutrinos (for the detailed calculations see [HI]), 
cr cff = (3tt/V2)Gf = 3.0 x 10~ 32 cm 2 is the effective 
cross-section obtained by integration of the Breit-Wigner 
cross-section over energy. 

The flux of i/ e -neutrinos produced in P7 cm b collisions 
is strongly suppressed and appears mostly due to oscilla- 
tions as has been already anticipated in (37| . According 
to recent calculations, Jp e (E) is approximately the same 
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for all flavours, namely it is about 1/6 of the all-flavour 
neutrino flux. 

For the application of the Glashow resonance to the 
Pop III bright phase, we have to first determine Zb- It 
can be found from the condition that the energy of the 
parent proton for the resonant neutrino at Zb must be 
above the critical energy E c /(1 + Zb), at which the energy 
losses for pion and e + e~-pair production are equal. Using 
E c = 6 x 10 19 eV one has 



( E c \ 1/2 
1 + Zb ={^) = 2L8 < 



(24) 



which coincides with our assumption Zb = 20. 

The attractive feature of this model is that it imposes 
milder constraints on E^ x . Indeed, the energy of the 
parent proton for resonant neutrino at z& is E p {zb) — 
20E a (1 + z b ) = 2.8 x 10 18 eV, which is considerably 
lower than E^ K discussed in section [TTJ Below we as- 
sume £ ac a c x = 1 x 10 19 eV. 

We calculate now the rate of the resonant events in 
IceCube with effective mass M ~ 1 x 10 9 t. According 
to Eq. (|23p . there is only one unknown quantity for the 
calculations, the flux Jp e (Eo), which we take as maxi- 
mal, namely corresponding to the cascade upper limit 
w™ = 5.8 x 1(T 7 cV/cm 3 . Using Eq. (JTHJl for the neu- 
trino flux with uj p (zb)/(l + Zb) 4 from Eq. (|2"0"j) . and as- 
suming equipartition of neutrino flavours with a fraction 
1/6 for u e , we obtain the rate of the resonant events in 
IceCube as 



0.1c 



N P a r 



12 ln(£™ ax / £p 



E 



(25) 



In more accurate calculations, where one takes into ac- 
count that decaying pions transfer only half of their en- 
ergy into e-m cascade, _E™ ax in Eq. ([25]) should be sub- 
stituted by ■J i?^i?™ ax , and the rate of events reaches 
v Ies = 1.3 yr _1 . In fact, this rate is not the maxi- 
mal one, because as already mentioned above, w™ ax = 
5.8 x 10 -7 eV/cm 3 was obtained for z = 0, while for 
Zb = 20 this limit is higher. Most probably each resonant 
event with tremendous energy release 6.3 x 10 15 eV in the 
form of nuclear and electromagnetic cascades will be de- 
tected in IceCube. Ten well identified events with equal 
energies during five years of observations is a good enough 
signature of Pop III bright phase model with Zb > 20. 



V. CONCLUSIONS 

There are theoretical and observational indications 
that the first star formation, the galaxy formation and 
AGN went through a phase of enhanced activity at high 
redshifts, z ~ 10 — 20 (e.g. low metallicity stars and the 
observation of about 20 quasars at z ~ 5 — 6.4 f2lj|). 

In this paper we concentrate upon a bright phase in 
the stellar evolution, in the form of Pop III stars. These 
astrophysical objects play an important role in at least 



two ways: 1) they enrich the universe with the metals 
that appear to be necessary for the formation of Pop 
II stars, thereby bridging the gap between the chemical 
composition of the universe as predicted by BBN and the 
one observed in nowadays universe. 2) Pop III stars may 
reionize the universe, at a redshift that WMAP observa- 
tions measure as z re ion = 11.0 ± 1.4 

Both the metal enrichment and the reionization require 
that Pop III stars are short-lived and therefore very mas- 
sive stars, typically with M* > 100M Q , which explode as 
supernovae. In the mass range 140M Q < M* < 260 Af© 
the star is fully disrupted due to the e + e~ pair-instability. 
In the mass range outside this interval a massive black 
hole and disc are produced, with ejection of gas in the 
form of jets. In both cases a SN explosion results in a 
shock where particle acceleration may take place. In the 
low-mass limit M ~ (15 — 40)M Q the energy release is 
W S n ~ 2 x 10 51 erg [33j], at large masses M > 100M Q 
it may exceed 10 53 erg. In case of a very massive black 
hole the energy output can be much higher due to the 
Blandford-Znajek effect. 

Here we concentrated on the possibility that diffusive 
particle acceleration may take place at the shocks that 
are generated when Pop III stars explode. The medium 
in which these SN explosions occur is expected to be 
made mainly of hydrogen (75 %), with a small contam- 
ination of helium and other light elements. Although 
a metal contamination must appear and increase once 
the bright phase starts, we expect that this represents 
a weak effect. Therefore, we considered protons as the 
main component of accelerated particles. 

We discussed at length three delicate aspects of the 
problem of acceleration: 1) the formation of collision- 
less shocks around primordial stars; 2) the generation of 
magnetic field; 3) particle scattering and acceleration. 

Pop III stars are the most natural sources of reioniza- 
tion of the universe due to their high temperature and lu- 
minosity fl6| . The reionization of the universe is a crucial 
ingredient of our discussion, since in the absence of such 
a phenomenon the collisionless shocks asociated with SN 
explosions would hardly form or would have very pecu- 
liar characteristics. These shocks are formed because of 
the mediation of electromagnetic instabilities such as the 
Weibel instability that was found to be effective even in 
the absence of a pre-existing magnetic field. The gen- 
eration of magnetic field is just a different aspect of the 
problem of shock formation, in that the presence of mag- 
netic field is the very reason why particle motion is slowed 
down, thereby fulfilling the Rankine-Hugoniot conditions 
at the shock surface. Magnetic field can also be formed 
downstream of the shock if the upstream plasma con- 
tains density inhomogeneities that induce shock corruga- 
tion and eventually lead to eddies that may considerably 
amplify a magnetic seed [40]. For relativistic shocks, am- 
plication of the magnetic field by a macroscopic turbulent 
dynamo triggered by the Kelvin-Helmholtz shear insta- 
bility has been investigated in |4lj . 

It is worth stressing that in these cases the magnetic 
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field is formed behind the shock surface and advected 
downstream. It is therefore not useful in terms of scat- 
tering the particles upstream of the shock, a necessary 
condition for particle acceleration. Magnetic field can be 
produced upstream of the shock through streaming in- 
stability excited by accelerated particles. This instabil- 
ity can be excited in a resonant or non-resonant way, the 
former being in general slower for fast shocks but more 
effective in scattering accelerated particles. We found 
that these processes can hardly allow particles to reach 
energies in excess of 10 19 eV for fast but non-relativistic 
shocks, although higher E max can be possibly obtained if 
numerous SNe explode in the same region. 

Higher maximum energies, in excess of 10 20 eV have 
been widely discussed and claimed to be achievable in 
the case of relativistic motion of the SN ejecta [13, [29| 
that are thought to occur in GRB afterglows (see however 
discussion in section ITT1) . 

UHE neutrinos in our model are produced in P7 cm b 
collisions. The energy of CMB photons is (1 + Zb) times 
higher than at z = 0, and the density of these photons 
is (1 + Zb) 3 times larger. We assume that a burst of 
UHE proton generation occurs at redshift Zb and con- 
sider two values Zb = 10 and Zb = 20 as benchmark cases. 
The generation spectrum is assumed to be cx E~ 2 with 
i?max = 1 X 10 20 eV or somewhat less. With these as- 
sumptions, the neutrino flux, given by Eq. Il~6"|) . is fully 
determined by the basic parameter represented by the 
energy density u p (zb) of UHE protons at epoch Zb, re- 
calculated at the present epoch ui p = uj p (zb)/(l + Zb) A 
(see Eq. (1161) ). In order to have the muon- neutrino 
flux detectable by IceCube, this parameter must be 
uj p > 1 x 10 -6 eV/cm 3 . The dominant number of de- 



tected events is confined within a limited energy interval 
centered at energy E c v = 7.5 x 10 15 (20/z 6 ) 2 eV (Eq. p^j) . 
At this energy the detected flux has a weak maximum. 
The predicted flux with the value of lu p given above is 
detectable by IceCube and respects the cascade upper 
limit on cosmogenic neutrinos and SN energetics for Pop 
III stars. 

A unique signature of the Pop III burst model may 
be provided by v e neutrinos interacting in a detector 
through the resonant reaction D e + e~ — > W~ — > hadrons. 
The resonant energy of the neutrino is given by Eq. (|22[) 
as Eo = 6.3 x 10 15 eV, and the rate of neutrino events 
in a detector is determined by one unknown quantity, 
the neutrino flux at energy E , J„ e (E ), see Eq. (|2"3"1) . 
We take this flux as to saturate the cascade upper limit 
corresponding to the cascade energy density w™f x = 
5.8 x 10 -7 eV/cm 3 , though in the model with fixed Zb this 
limit is expected to be higher. The redshift of the burst 
Zb, or more precisely its lower limit, can be found from 
the condition that the energy E p (zb) of a proton, parent 
of the resonant neutrino, must have energy higher than 
E c (zb)i & t which pion energy losses at epoch Zb exceed 
those for e + e~ production. This condition determines 
the burst redshift as Zb — 20.8, see Eq. (j24|) . The resonant 
events can be observed in IceCube through the Cercnkov 
light from nuclear and e-m cascades with energy Eq. In 
IceCube the predicted frequency of resonant events for 
the flux at the cascade upper limit is rather low, about 
10 events in 5 years, but with a tremendous energy de- 
posit, the same for all events. This may provide us with a 
reasonable signature of the model discussed here. These 
events are accompanied by muons with energies of order 
0.5££ ~4x 10 15 eV. 
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